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A model describing the three-dimensional folding of the triangular lattice on the face-centered 
cubic lattice is generalized allowing the presence of defects, which are related to cuts in the two- 
dimensional network. The model can be expressed in terms of Ising-like variables with nearest- 
neighbor and plaquette interactions in the hexagonal lattice; its phase diagram is determined by 
means of the Cluster Variation Method. The results found by varying the curvature and defect 
energy show that the introduction of defects turns the first-order crumpling transitions of the model 
without defects into continuous transitions. New phases also appear by decreasing the energy cost 
of defects and the behavior of their densities has been analyzed. 
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modynamics of surfaces and interfaces) 
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I. INTRODUCTION 

The behavior of fluctuating membranes and surfaces is 
relevant for many physical and biological systems, from 
gauge theories and strings to vesicles and cellular mem- 
branes; see Refs. [TJ [2] for reviews. A class is given by 
polymerized or crystalline membranes [3H5] consisting of 
two-dimensional networks of molecules with fixed con- 
nectivity. Examples are the spectrin network in red- 
blood cells [5] or the graphene [7]. Theoretical argu- 
ments suggest the relevance of curvature energy terms 
for the macroscopic behavior of fluctuating surfaces pQ. 
For polymerized membranes, at variance with fluid mem- 
branes [8J, a flat phase with long-range order in the ori- 
entation of the normals to the surface is expected to be 
stable at high bending rigidity [9]. 

The prediction of a crumpling transition for phan- 
tom polymerized membranes (where self-avoidance is 
not taken into account) is confirmed on the basis of 
various analytical results 1 1 Ofi 1 31 and numerical simula- 
tions [3J [TS] on continuous models. The character 
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of the transition is controversial. Among the most re- 
cent results, the simulations of Refs. suggest a first- 
order behavior while the non-perturbative renormaliza- 
tion group calculations of 17] predict for a phantom 
membrane embedded in the three-dimensional space a 
continuous transition. 

A different approach has been to study discrete models 
for polymerized membranes whose nodes are constrained 
to occupy positions corresponding to the sites of a given 
lattice. Solutions of discrete models, as in other domains 
of statistical physics, can represent a reference for the 
behavior of fluctuating membranes of a given class. It 
is well-known that Ising-like models [TSJ [TU] naturally 
admit an interpretation in terms of a surface gas with 
a variable number of components. More difficult is the 
problem of representing a single membrane in terms of 
discrete variables. In Ref. 20J Bowick and co-workers 
have been able to model the folding of a triangular net- 
work in the face-centered cubic (fee) lattice in terms of lo- 
cal discrete variables, subject to local constraints. Since 
the distance between the nodes of the network is fixed, in 
this model one can consider as frozen the "phonon" de- 
grees of freedom of the membrane while only the bending 
modes are taken into account |21j . A simpler version of 
the folding problem of [20] has been studied in [22H25] 
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corresponding to the case of a two-dimensional embed- 
ding space. Here the normals to the triangles of the net- 
work can point only "up" or "down" i n some direction. 
A similar model concerning the folding properties of a 
square lattice along the main axis and the diagonals has 
been studied in j^Hl HZ] . Finally, the three-dimensional 
folding problems of the triangular lattice with quenched 
random bending rigidity and spontaneous curvature have 
been also respectively studied in [28] and [25] , 

The phase diagram of the model of [5D] has been first 
studied in [3D] and [3T] by the cluster variation method 
(CVM) [32] ■ A sequence of transitions from the flat to the 
piled-up phase through the partially folded octahedral 
and tetrahedral phases was found by varying the bending 
rigidity from oo to — oo. The density-matrix renormal- 
ization group (DMRG) calculations of [33J [34] have con- 
firmed the occurrence of these transitions. Both CVM 
and DMRG have predicted first-order flat-octahedral 
and octahedral-tetrahedral transitions while the charac- 
ter of the third transition is controversial. CVM sug- 
gested a continuous tetrahedral piled-up phase transition 
at variance with the weak first-order behavior predicted 
by DMRG. 

In order to further analyze the character of the crum- 
pling transitions of the triangular network in the fee lat- 
tice, a modified folding rule allowing defects in the net- 
work has been introduced in [35] . Transfer-matrix cal- 
culations of [35] confirmed the order of the transitions 
given by DMRG. 

In the present work we observe that the model of fold- 
ing with defects is interesting by itself. Relaxing the 
local constraints on the Ising variables in the model of 
Ref. [20] corresponds to accepting folding configurations 
with bending and/or meeting points in the folding lines, 
which can be obtained by allowing cuts between adja- 
cent triangles. Therefore the model with progressively 
relaxed constraints can describe polymerized membranes 
with an increasing number of defects in the connectivity 
rules. Here we will introduce such a model and study its 
equilibrium properties by means of the CVM. 

The paper is organized as follows. In the next sec- 
tion we will introduce the model and describe the CVM 
approximation scheme used for studying it. Our results 
for the phase diagram and the equilibrium behavior of 
the most important quantities are reported in Sect. III. 
Conclusions will complete the paper. 



II. THE MODEL AND THE METHOD 

The model for the folding of the triangular lattice in 
the three-dimensional fee lattice has been introduced in 
[2pJ[3pJ[3l]. 

In this model, sites of the triangular lattice are mapped 
onto sites of the fee lattice, with the condition that two 
sites that are nearest-neighbors (NNs) in the triangular- 
lattice must remain NNs in the fee lattice. Let us consider 
two adjacent plaquettes of the triangular lattice, and call 
9 the angle formed by their normal vectors (defined in 
such a way that 9 always vanishes in the planar configu- 
ration) . It is easy to see that upon mapping onto the fee 
lattice, the two plaquettes will be in one of four relative 
orientations, namely: (i) no fold, with 8 = 0,cos# = 1; 
(ii) octahedral fold, with cos6> = 1/3 (the two plaquettes 
belong to the same octahedron in the fee lattice); (iii) 
tetrahedral fold, with cos 8 = — 1/3 (they belong to the 
same tetrahedron in the fee lattice); (iv) complete fold, 
with cos 9 = — 1 (they are on top of each other). A fold 
between two adjacent plaquettes has an energy cost, due 
to curvature, given by — K cos 9. 

The model allows various representations |20) . a par- 
ticularly convenient one for our purposes being defined 
in terms of 2 sets of Ising variables. These variables are 
usually denoted by Oi and Zi, the pair (<7j, Zi) being asso- 
ciated to a plaquette of the triangular lattice or, equiva- 
lently, to a site i of the dual hexagonal lattice A. Sites of 
A correspond to plaquette centers in the triangular lat- 
tice, while edges of A are perpendicular to the edges of 
the triangular lattice. In terms of these Ising variables, 
the fold between two adjacent plaquettes 1 and 2 (corre- 
sponding to two NN sites in A) is specified by the values 
(7i<72 and Z\Z2 according to Tab. [TJ and the angle 8 is 
given by 



In the language of magnetism, we can say that different 
folds correspond to different types of domain walls for 
our Ising variables. 

However, not all configurations of our Ising variables 
are allowed. It has been shown in [20] that the variables 
Zi and <Tj have to satisfy two constraints, or folding rules, 
in order to describe a proper folding configuration over 
the fee lattice. Given an hexagon e in the set E of all 
the elementary plaquettes in A, let 1,...,6 be its six 
sites ordered counterclockwise, a e — {<Ji, i = 1, . . . 6} and 
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Type of fold 


cos 6 


0"1<72 


Z\Zl 


no fold 


+ 1 


+ 1 


+ 1 


octahedral 


+1/3 


-1 


-1 


tetrahedral 


-1/3 


+ 1 


-1 


complete fold 


-1 


-1 


+ 1 



TABLE I. Type of folds and their representation in terms of 
Ising variables. 
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z e = {zi, i = 1, . . . 6}. Set 

6 

and 

M c (cr e , z e ) 
where z^ — Z\ and 
A i)C (cr e ) = 



(2.2) 



(2.3) 



i=l 



1 if J2j=i <J j — c mod 3 
otherwise 



with i = 1, 
read 



(2.4) 

6 and c = 1.2. The folding rules then 



L(tj e ) = mod 3 
M c (a e ,z e ) = mod 2, 



(2.5) 

1,2,3. (2.6) 



Notice that M 1 (a e , z e )+M 2 ((Je, z e )+M 3 (a e , z e ) is an even 
number, and therefore the number of violated constraints 
in 2.6 can be only or 2: as a consequence, it is sufficient 
to impose only 2 of the 3 conditions in |2.6| 

We now introduce defects in the configurations of 
the triangular lattice. A violation of the folding rules 



Eqs.(2.5), (2.6) at a hexagon represents a defect at that 



hexagon. 



The folding rules Eqs.(2.5) and (2.6) ensure that after 



a complete turn around any set of six triangles with a 
common vertex the triangles occupy the same absolute 
position [2D] . Clearly, defects alter the local connectivity 
of the original triangular lattice and can be associated 
to cuts in the lattice. More specifically, defects can be 
associated to certain kinds of bending and/or meeting 
points of the various domain walls for the (cr, z) variables, 
which are not allowed in the original folding model and 
therefore do not appear in the set of allowed vertices in 
Fig. 10 of [20] • 

In Fig.JT] for example, we have a hexagon configuration 
containing only a domain wall between (cr, z) = (+1, +1) 



and (— 1, +1) (we call it a complete-fold domain wall), 
which shows a bend at the center of the hexagon. This 



corresponds to a violation of constraint (2.5 1 only, as one 
can easily check. This configuration cannot appear in 
a fold of the triangular lattice, unless one allows cuts. 
It is however immediately evident that, at least locally, 
there are three possible ways to cut the hexagon. One 
can select one of the three edges of triangle number 6, 
cut the remaining two edges, and then fold the triangle 
along the selected edge. If this were the only defect in the 
whole configuration, this argument would apply globally 
as well. Therefore, already on the basis of this example, 
it is clear that a unique reconstruction of the surface is 
not feasible in the presence of defects. 



The analysis of a similar configuration with ((Ti,Zi) = 
(+1, +1), i = 1, . . . , 5 and (a§, z 6 ) = (—1, —1), represent- 
ing the bending of an octahedral-fold domain wall shows 
that in this case both constraints (2.5) and (2.6) are vi- 
olated. On the other hand, in the case of bending of a 
tctrahedral-fold domain wall between (cr, z) = (+!,+!) 



and (+1,-1), only (2.6) would be violated. On the ba- 
sis of the previous examples one might be tempted to 
consider a specific geometrical interpretation for the sep- 
arate violation of each of the two folding rules. This 
other example, however, shows that this interpretation 
would be problematic. In Fig. [2j we see a configuration 
where an octahedral-fold domain wall meets (ending at 
the meeting point) another domain wall, which at the 
meeting point changes its character from tetrahedral to 
complete: in this case both constraints (2.5) and (2.6) 
are violated. If, however, the complete and octahedral 
domain walls are swapped, only (2.5) is violated. One 
can conclude from the above examples that the two fold- 
ing rules can not be directly associated to definite classes 
of defects. 



Due to the above observation, we do not find reasons 
to attribute different energy costs to different kinds of 
defects and will consider a model where all defects are 
weighted in the same way [36] . A single energy parame- 
ter A will be coupled to the number of defects, that is the 
number of hexagons at which (one or more) constraints 



are violated. Taking also into account Eq. (2.1 ) for writ- 
ing the curvature energy as in 



we are led to consider 
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the Hamiltonian (energies are given in units of k^T) 



K 

H = — - 22^^(1 + 2z i z j ) 



eG_E 



{L(<r e )=0 mod 3} 
3 

W I{A/ c (ct c ,z c )=0 mod 2} 
c=l 



(2.7) 



where I{ CO ndition} is equal to one if the condition is sat- 
isfied and to zero otherwise. Note that in the above 
formula the first sum is extended to the NN pairs and 
the second to the hexagons of the lattice. Moreover, the 
product over the three color indices in the last sum is re- 



dundant: according to our discussion of constraints (2.6 1, 
it is sufficient to impose only 2 of the 3 constraints. 

Finally, we observe that it is not possible to express 
the energy cost of defects in terms of local weights for the 
length of the cuts needed to obtain a surface realization of 
a given (er, z) configuration. This can be shown by the fol- 
lowing example. Consider a set of 3 defects correspond- 
ing to 60° bends of complete-fold domain walls as those 
shown in Fig. [T] These defects can be used to construct 
configurations where a complete-fold domain wall can 
form equilateral triangles of any size, (er, z) = (+1,+1) 
outside the triangle, and (<r,z) = (— 1, +1) inside it. 
Hence, given the same set of 3 defects, one can construct 
configurations corresponding to cuts of any size. 

As it was done for the two-dimensional model dis- 
cussed in [24] and for the defect free three-dimensional 
case studied in [3D], the phase diagram will be inves- 
tigated by means of the hexagon approximation of the 
CVM. For the sake of self-containedness we briefly recall 
the main features of our approach. 

We recall that E is the collection of hexagons on 
the lattice. We denote by M the collections of all the 
hexagons and all their sub-clusters (site subsets). To 
each cluster a G M a probability distribution p a (a a , z a ) 
is associated, where a a = {(Xi,i G a} and z a = {zi,i G 
a}. In this CVM approximation the free energy func- 
tional is given by 



F * = E E E ^ a 

a : Z& ) Pa (&at : %a ) 



aeM T a ,z a 



2J a a ^2 Pa{(Ja,Z a )\0g 

Pa [Gat j %ct ) •> 



(2- 



where the coefficients a a are such that for each a G M 

E a p = 1 ( 2 - 9 ) 

/3GM:/3Da 

and the Hamiltonian terms H a are defined as follows: for 
a cluster made of two neighboring sites we have 



H (i,o) = 



K 



<7iOj{\ + 2z i z J ), 



for a cluster made of a hexagon ee£ 



H,. 



-XI 



{L(ct c )=0 mod 3} \\ J1 {A/ c (cr c ,z e )=0 mod 2} 
c=l 



and H a = otherwise. By using Eq. (2.9 1 it is not diffi- 



cult to prove that the sole non-vanishing coefficients a a 
are those associated to hexagons, NN pairs, and single 
sites; coefficients associated to five-site, four-site, three- 
site, and not neighboring two-site clusters are all equal 
to zero. In particular we have that 



a e = 1, a 



-1, and a,i = 1 



for each hexagon e G E, each nearest-neighbor pair 
and each site i G A. 

We now exploit the translational invariance of the sys- 
tem, which implies that all the probability distributions 
associated to a particular sub-family of clusters in M are 
equal. Thus, we denote by pe(ai, . . . , eg, Zi, . . . , Zq) the 
probability distribution associated to hexagons and by 
P2(ci, 02, Z\, z%) that associated to NN pairs. Site clus- 
ters form two sub-families: indeed, a hexagonal lattice 
is a bipartite lattice, made of two inter-penetrating (tri- 
angular) sub-lattices a and b, such that all the NNs of 
a site in a belong to b and vice-versa. Since we expect 
the symmetry between the two sub-lattices to be broken 
in some thermodynamic phases of the model, it is im- 
portant to distinguish the site probability distributions 
corresponding to a and b. For the same reason, in the ar- 
gument of p 2 , the first and the third entries (cti and Zi) 
refer to sub-lattice a, while the second and the fourth 
ones ((72 and z 2 ) refer to sub-lattice b, and they cannot 
be interchanged. 

With the above definitions and notations, and observ- 
ing that the number of hexagons and NN pairs in A are 
respectively N 6 = |A|/2 and N 2 = 3|A|/2, from Eq. Kty 



we obtain the following expression for the CVM free en- 
ergy density functional / = Fjvf/|A|: 



aeM <r a .z a 
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1 1 

/ = --KTr 2 [aia 2 {l + 2ziZ 2 )p 2 ] - ~Tr 6 \Xl{L((T e )=0 mod 3} JJ I{M c (o- e ,z e )=0 mod 2})pe 

c=l 



(2.10) 



13 1 1 

-Tr 6 [p 6 lnp 6 ] - -Tr 2 [p 2 ln/j 2 ] + -Tri^i^ lnpi >a ] + -Th In p lj6 ] + ^(Tr 6/ 9 6 - 1) 



where we have introduced the notation 



T l ."2 

si ,^2 



and i/ is a Lagrange multiplier which ensures the normal- 



r 



ization of p&. The other probability distributions do not 
need a normalization constraint since they can be writ- 
ten as partial traces of probability distributions of larger 
clusters. More precisely, 



g E lP$( a l> a 2i a 3> a 4> a 5, f6i ■ ■ ■ ) + P6{<^3, °2, CTi, cr 4 ,cr 5 , CT 6 , . . . ) + p 6 (a 3} a 4: , (Ti, cr 2 , CT 5 ,cr 6 ,...) ( 2 -H) 

CT 3 i ■ ■ ■ ' CT 6 

+P6(c r 3 I cr 4 ,a 5 ,a 2 ,a 1 ,a 6 , . . .) + p 6 (a 3 ,a4,a 5 ,a 6 ,a 1 ,a 2 , . . .) + p 6 (a 1 ,a 4: ,a 5 ,a 6 ,a 3 ,a 2 , . . .)} 



where the z variables in the argument of pe appear in the 
same order as a variables, and 

Pl,a(<ri,zi) = E P2(o-i,ct 2 ,zi,z 2 ), (2.12) 
/Oi,b(cr 2 ,z 2 ) = E P2(o'i,cr 2 ,Zi,z 2 ). (2.13) 

fl,Zl 

With the above definitions the CVM free energy den- 
sity functional / can be regarded as a function of p^ 



I 

only. The minimization must be performed numerically, 
and this can be easily done by standard iterative meth- 
ods as in [24] (see [32] for a survey of such algorithms). 
The simplest possibility is to write stationarity equa- 
tions by taking derivatives of / with respect to an ele- 
ment of Pq{(Ti, . . . , cr 6 , z±, . . . , Zq), for some generic choice 
<7i, . . . , eg j z ii ■ ■ ■ i z 6 — ±1 of the spin variables and let- 
ting (T7 = (7i and Z7 = z\ . After some algebra, we get 



p§(p\,. . . ,cr 6 ,zi, . . . ,z 6 ) 

6 



= exp{(K/6) E OiOi+iC 1 + 2ziZi+i) - 2u + AI {L(CTs)=0 mod 3} J| '. 



■{M a [<7 e ,z e )=0 mod 2} 



} 



(2-14) 

x [p 2 {a 1 ,a 2 , zi,z 2 )p 2 (a 3 , a 2 , z 3l z 2 )p 2 (a 3l cr 4 , z 3 , z 4 )p 2 (a 5 , cr 4 , z 5 , z 4 )p 2 {a 5 ,a 6 , z 5 , z 6 )p 2 (<ri, <J 6 , z\, z 6 )} 1/2 

X [pi,a(<7l, 2l)pi,b((7 2 , Z 2 )pi : a((7 3 , Z 3 )Pl,f)(o-4, 2*0pi,a ("5 , Z 5 )pl,b (^6 , Z 6 )]~ 1/3 



that can be solved numerically with an iterative ap- 
proach. 



III. RESULTS 



We now present our results for the phase diagram of 



the folding model Eq. (2.7|, obtained by finding the sta- 



ble (lowest free energy) solutions of Eq. (2.14) at varying 



curvature energy and defect cost. Following |31j we in- 
troduce the order parameters O = (erf), T = (zjcrf*) and 
P = (erf*). O and T were named octahedral and tetrahe- 
dral order parameters, respectively, and indeed they are 
the only order parameters which become non-zero in the 



G 



Phase 


M 


p 


N 


o 


T 


flat 














p-octahedral 
















tetrahedral 














7^0 


piled-up 












7^0 


f-octahedral 
















disordered 


















TABLE II. Phases. In the first column the name of the phase 
is reported. In the columns from the second to the sixth it 
is indicated which, among the order parameters M = (en), 
P = (zi), O = (erf), N = {ana), and T = {afm), differ from 
zero. 

corresponding phases (observe, as explained later, that 
the octahedral phase of the model without defects will be 
called in this paper p-octahedral). P was named planar 
order parameter, since it becomes non-zero in the pla- 
nar (here called flat) phase. In order to characterize and 
distinguish all the phases we obtain here, it is convenient 
to define two additional order parameters, M = (crj) and 
N = (zitJi). M = 1 corresponds to configurations with 
tetrahedral folds or no folds only, while N = 1 corre- 
sponds to configurations with octahedral folds or no folds 
only. M is non-zero only in the flat phase, while N is 
non-zero in the flat and in the f-octahedral phase, as 
described later. The identification of the various phases 
of the model in terms of the above order parameters is 
summarized in Tab. ITT1 

Fig. [3] summarizes our results for the phase diagram in 
the plane K-X. We shall describe it by first considering 
the large A limit, where defects are absent and one re- 
covers previous results [30] [31] . We shall then consider 
smaller values of A in order to see how the various phases 
and phase transitions are modified by the introduction of 
defects. 

At large enough A, folding rules are practically never 
violated and the equilibrium thermodynamics of the 
model becomes independent of A. Indeed, in the upper 
portion of the phase diagram Fig. [3j the phase transition 
lines are practically vertical and do not change their na- 
ture anymore. We can easily check that the results of 
the defect-free case [30l EI] are recovered in this limit. 
For instance, at A = 20 and vanishing K we obtain the 
entropy (per site) S = lng, where q — 1.42805, equal to 
the value obtained by the constrained CVM approach of 
[30j and in very good agreement with the transfer ma- 
trix estimate q = 1.43(1) of [20]. Furthermore we obtain 



O = 0.87456, M = P = N = T = 0, indicating a 
marked preference of the triangular lattice for wrapping 
on an octahedron at zero curvature cost. As already no- 
ticed in [31] . here most folds between adjacent plaque- 
ttes are octahedral or complete. The sequence of folding 
transitions observed in the defect-free limit J3U1 EI] is 
reproduced here already at A > 4.91. Choosing A = 10 
as an example, we find a flat phase at large K (in the 
limit A — > 00 this phase is perfectly flat in this approx- 
imation), then at K = 0.1856 a first-order transition 
occurs between the flat and the octahedral phases, char- 
acterized by M = 0, P = 0, N = 0, O = 0.8247, and 
T = 0. This phase, as shown in [31], is also characterized 
by a relative abundance of complete folds with respect to 
other folds. In order to make a difference with another 
phase appearing in the phase diagram (see the following) 
that is characterized by a relevant presence of octahedral 
folds, it will be called p(iled-up)-octahedral phase. Upon 
further decreasing K we find another first-order transi- 
tion at K = —0.2940 between the p-octahedral phase 
with M = 0, P = 0, N = 0, O = 0.5816, and T = 
and the tetrahedral one with M = 0,P = 0,N = 0, 
O = 0, and T = 0.7466, Finally, at K = -0.8395 there 
is a continuous transition from the tetrahedral phase to 
the piled-up phase with the order parameters charac- 
terized by the continuous vanishing of P and O, while 
T = 0.9993. This sequence of phases and phase transi- 
tions agrees in nature with that found in the defect-free 
limit [30] EI] 1 providing a confirmation of the validity of 
the present approach. 

Moving to lower values of A, for A > 4.91 the behavior 
described above remains qualitatively the same as in the 
case without defects, as it can be seen in Fig. [3j An 
illustration, for A = 6, in terms of the order parameters as 
functions of K, is reported in the right column of Fig. [4] 

By further decreasing the value of A a new phase, 
characterized by the vanishing of all the order param- 
eters, appears in the phase diagram. This phase, which 
we shall call disordered, is bounded by the tetrahe- 
dral and the p-octahedral phases at negative K and by 
the p-octahedral, the flat and the f-octahedral (see be- 
low) phases at positive K . The p-octahedral-disordered 
and tetrahedral-disordered phase transitions are contin- 
uous and the corresponding critical lines meet the p- 
octahedral-tetrahedral first-order transition line at a bi- 
critical point at K = -0.298, A = 4.90 (see Fig. [3f . At 
positive K, the flat-disordered transition is first-order, 
like the flat-p-octahedral one, and the p-octahedral- 
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disordered critical line meets them at a critical end-point 
at K = 0.208 and A = 3.84. At zero curvature, the 
p-octahedral-disordcrcd transition occurs at A = 3.478. 
The behavior of the order parameters as a function of K 
at A = 3.6 is reported in the central column of Fig. [4] 

Let us now proceed by describing our phase diagram at 
positive K and small A. In this region another new phase, 
to be denoted by f(lat)-octahedral, is found between the 
disordered and the flat phase. The f-octahedral phase is 
characterized by N being the only non-vanishing order 
parameter, while in terms of folds between adjacent pla- 
quettes it exhibits a mixture of no-folds and octahedral 
folds. This can be seen in Fig. [5] where the average pro- 
portions of the four types of folds is shown as a function 
of K, for three values of A, A = 0, 3.6 and 6. 

The appearance of the f-octahedral phase can be bet- 
ter understood by discussing the behavior of the model 
at A = 0. As observed in [31], on this axis the model re- 
duces to the Ashkin-Teller model [37], with parameters 
corresponding to the trivial case of two independent Ising 
models. Starting from low K, two continuous, symme- 
try breaking transitions are found. The first transition, 
separating the disordered phase from the f-octahedral 
phase, occurs at K — 3K c /2, where K c is the critical 
Ising coupling on the hexagonal lattice, as discussed in 
[31j . Here the order parameter N becomes different from 
zero. The second transition, separating the f-octahedral 
phase from the flat phase, occurs at K = 3K C [31]. Here 
the symmetries z — z and a — > — a are separately bro- 
ken, and the order parameters M and P also become 
different from zero. 

We find, in our approximation, that the two transi- 
tions occur at K c \ = 0.9321 and at K C 2 = 1.8642. These 
values correspond to the CVM estimate of the Ising crit- 
ical point K c — 0.6214 on the hexagonal lattice (in the 
hexagonal plaquette approximation) 24 . At positive A 
the transitions separating the f-octahedral and the dis- 
ordered phases and the f-octahedral and the flat phases 
remain continuous. The corresponding lines meet at a 
bicritical point at K — 0.388, A = 2.39 where the first- 
order disordered-flat transition appears. 

At negative K the tetrahedral-piled-up phase transi- 
tion remains continuous for all values of A. The transition 
line intersects the horizontal axis at K = ~K C 2 where 
the symmetries z — > — z and a s t — > — cr st are separately 
broken. Moreover, at A = 0, the disordered-tetrahcdral 
continuous transition is found at K = —K c \- These criti- 
cal values correspond to the anti-ferromagnetic images of 



the Ising transitions occurring in the equivalent Ashkin 
Teller model at positive K . The behavior of the order 
parameters as a function of K at A = is reported in the 
left column of Fig. [4] 

The above results show that the introduction of de- 
fects turns first-order transition into continuous one, as 
in the planar folding case [21] ■ Here however the phase 
behavior is richer and we can observe an additional effect 
due to the presence of defects: the p-octahedral phase 
disappears, to be replaced by the disordered phase and 
the f-octahedral phase, which also exhibits a significant 
fraction of octahedral folds. 

It is also worth taking a look at the number of defects 
appearing in our triangular lattice as a function of K 
and A. To be more precise we define the fraction pl of 
hexagons at which folding rule Eq. (|2.5|) is violated 



PL = Tr 6i 6 (l - I{L( C r e ) =0 mod 3}) > (3.15) 

the fractions pm c of hexagons at which folding rules Eq. 



(2.6) are violated 



PM C — Tr 6 p 6 (l - I{M c (ve,Ze)=0 mod 2}) C — 1, 2 

(3.16) 

and the fraction p of hexagons at which at least one fold- 
ing rule is violated, 

V = Tr 6 pe(l - H{M o (o- e ,z e )=0 mod 2} x 

-A, \ ( 3 - 17 ) 

X [ [ I{M o (>e,Ze)=0 mod 2} )• 
c=l 

These quantities are plotted as functions of K in Fig.[6]for 
A = 6 (that is, close to the defect-free limit), in Fig.[7]for 
A = 3.6, and in Fig.[8]for A = 0. We see that the fraction 
of defects has a maximum close to K = and decreases 
as \K\ increases. Moreover, the fraction of defects is a de- 
creasing function of A. Considering the various phases we 
see that the phases exhibiting less fluctuations, namely 
the flat and piled-up phases, are almost defect-free, even 
at small A, while the other phases are more prone to de- 
fects. 

Notice also that pm x and pm 2 differ in the flat, p- 
octahedral and piled-up phases. This difference between 
Pm ± and p^ is a consequence of the breaking of the 
global inversion symmetry <Ji — ¥ —Ui,Zi — z,-,Vi. This 



transformation maps configurations violating (2.6) for 
c = 1 to configurations violating (2.6) for c = 2, and 



viceversa. The corresponding symmetry is preserved only 
in the disordered, f-octahedral and tetrahedral phases, 
where pm x = Pm 2 as a consequence. 
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IV. CONCLUSIONS 

We have generalized a model for the folding transitions 
of a triangular lattice in a three dimensional space, dis- 
cretized as a fee lattice, by allowing defects corresponding 
to cuts in the triangular lattice, and weighing them by 
a suitable energy cost. We have studied the model in 
a six-point approximation of the CVM. In the limit of 
the energy cost of a defect going to +00, we recover pre- 
vious results for the defect-free model [301 GH] ■ O n the 
other hand, when this energy cost is sufficiently small, 
we find that first-order transitions are turned into con- 
tinuous one, and that the octahedral phase found in the 



defect-free limit, which was characterized by octahedral 
and complete folds between adjacent plaquettes, is re- 
placed by a fully disordered phase and another phase 
characterized by octahedral folds and no-folds, which 
is related to the intermediate temperature phase of the 
Ashkin-Teller model. The model has a rich phase dia- 
gram with several multicritical points, namely, two bicrit- 
ical points and a critical end-point. We have also shown 
that defects are more likely to occur in phases exhibit- 
ing larger fluctuations, while the flat and piled-up phases 
are almost defect free, and their concentration typically 
decreases as the absolute value of the curvature energy 
increases. 
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2 4 = + l 
.0"4 = + l- 



/z 5 = + l\ 
OS=+l 

C6=— 1 . 



23=- 

<T2= I 1 . 

. 22 = 



2^ 



(7l= + l 



FIG. 1. On the left: a hexagon configuration with a complete fold domain wall represented by the thick line between (<j,z) = 
(+1.+1) and (— 1,+1). On the right: the three possible surface folded configurations corresponding to the hexagon spin 
configuration on the left, obtainable by cuts in the original triangular network as described in the main text. 
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FIG. 2. Two different hexagon configurations where a complete fold, a tetrahedral and an octahedral domain wall (respectively 
represented by lines of decreasing thickness) meet at the center of the hexagon. The octahedral and the complete fold domain 
walls are swapped in the two configurations. 




FIG. 3. Phase diagram of the model Eq. (2.7 1, where P, T, PO, D, FO, and F stand respectively for piled-up, tetrahedral, 
piled-octahedral, disordered, flat-octahedral, and flat. The solid circles and the solid box denote, respectively, the bicritical 
points and the critical end-point. Solid and dashed lines denote, respectively, first-order and continuous transitions. 
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FIG. 4. Order parameters vs. K for different values of the constraint parameter A. From the left to the right A takes the values 
0, 3.6, and 6. The thin vertical lines indicate points where a phase transition occurs. From the left to the right the phases 
are piled-up, tetrahedral, disordered, flat-octahedral, and flat at A = 0, piled-up, tetrahedral, disordered, piled-octahedral, 
disordered, and flat at A = 3.6, piled-up, tetrahedral, piled-octahedral, and flat at A = 6. 
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FIG. 5. Different fold fractions. From the left to the right A takes the values 0, 3.6, and 6. The thin vertical lines indicate points 
where a phase transition occurs. From the left to the right the phases are piled-up, tetrahedral, disordered, flat-octahedral, 
and flat at A = 0; piled-up, tetrahedral, disordered, piled-octahedral, disordered, and flat at A — 3.6; piled-up, tetrahedral, 
piled-octahedral, and flat at A = 6. The dotted lines, in increasing thickness order, denote respectively no fold, octahedral fold, 
tetrahedral fold, and complete fold. 
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FIG. 6. Fraction of defects as a function of K for A = 6. The dotted lines, in increasing thickness order, denote respectively 
Pl, PMi, Pm 2 , and p. The thin vertical lines indicate points where a phase transition occurs. From the left to the right the 
phases are piled-up, tetrahedral, piled-octahedral, and flat. 
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